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Abstract. The 7 s X-ray pulsar IE 2259-1-586 and the su- 
pernova remnant (SNR) G109.1-1.0 (CTB 109) were ob- 
served by BeppoSAX in 1996 November. The pulse period 
of 6.978914 ± 0.000006 s implies that IE 2259+586 con- 
tinues its near constant spin-down trend. The 0.5-10 keV 
pulse shape is characterized by a double peaked profile, 
with the amplitude of the second peak ~50% of that of 
the main peak. The pulse profile does not exhibit any 
strong energy dependence. We confirm the ASCA discov- 
ery of an additional low-energy spectral component from 
lE2259-f586. This can best be modeled as a 0.44 keV 
blackbody, but we cannot exclude that some, or all, of 
this emission arises from the part of the SNR that lies 
within the pulsar's extraction region. 

The spectrum of G109. 1 — 1.0 is well fit with a non- 
equilibrium ionization plasma model with a best-fit tem- 
perature of 0.95 keV. The derived mass for the X-ray emit- 
ting plasma ('^15-20Mq) and its near cosmic abundances 
imply that the X-ray emission comes mainly from mildly 
enriched, swept-up circumstellar material. The spectrum 
is strongly out of equilibrium with an ionization age of 
only 3000 yr. This age is in good agreement with that 
derived from hydrodynamic simulations of the SNR using 
the above X-ray temperature. 
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1. Introduction 

IE 2259-1-586 is an unusual pulsar with a pulsation pe- 
riod of 7 s and a very steep spectrum (photon index, a, 
~4). The source has exhibited a near constant spin-down 
trend at ~5xl0~^^ ss~^ since its discovery. IE 2259+586 
is located in, or near, the supernova remnant (SNR) 
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G109.1-1.0 which is also known as CTB 109. The SNR 
consists of an approximately hemispherical shell of X-ray 
and radio emission and a jet-like lobe located mid- way be- 
tween the shell and the pulsar (Gregory & Fahlman 1980; 
Hughes et al. 1981; Morini et al. 1988; Hurford & Fesen 
1995). There may be a compact (2-3' radius) synchrotron 
nebula surrounding the pulsar (Rho & Petre 1997). The 
location and shape of the X-ray lobe led to speculation 
that it is a "jet" of material connecting the pulsar and 
the SNR (Gregory & Fahlman 1983). However, its ther- 
mal spectrum and detailed morphology do not support 
this view (Hurford & Fesen 1995; Rho & Petre 1997). 

IE 2259+586 together with other sources such as 
IE 1048.1-593 and 4U 0142+614 belong to a small group 
of pulsars with similar spin periods (around 6 s) and prop- 
erties that clearly distinguish them from the "classical" 
pulsars in high-mass X-ray binaries (Mereghetti & Stella 
1995). Koyama et al. (1987) argue that the slow spin-down 
rate implies that IE 2259+586 contains an accreting neu- 
tron star rotating at close to its equilibrium period with a 
magnetic field of ~5 x 10^^ G. However, there is no direct 
evidence that the pulsar is in a binary system. In par- 
ticular, no optical or radio counterparts have been found 
(Coe & Jones 1992; Coe et al. 1994) and no Doppler shifts 
in pulse period have been detected (Koyama et al. 1989; 
Hanson et al. 1988; Morini et al. 1988). 

The lack of an optical counterpart and orbital Doppler 
shifts argue against a binary models for IE 2259+586, un- 
less the companion has an extremely low mass. Models 
based on single stars have been proposed, such as spin- 
down of a white dwarf (Paczynski 1990; Usov 1994), loss 
of magnetic energy of a strongly magnetized neutron star 
(Thompson & Duncan 1993) and a neutron star accreting 
from a circumstellar disk (Corbet et al. 1995; van Paradijs 
et al. 1995). The detection by Mereghetti (1995) of a 
large change in the spin-down rate of the similar system 
IE 1048.1— 5937, as well as the long term pulse frequency 
fluctuations in IE 2259+586 (Baykal & Swank 1996) sup- 
port an accretion hypothesis. 
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We present a detailed study of the X-ray emission 
from IE 2259+586 and G109. 1 — 1.0 using the imaging ca- 
pabilities of BeppoSAX to separate the contributions of 
the individual sources. This topic is particularly interest- 
ing because of the disparate spectral results reported for 
IE 2259+586 as well as of the uncertain origin of the X-ray 
lobe in G109. 1-1.0. 

2. Observations 

Results from the the Low-Energy Concentrator Spectrom- 
eter (LEGS; 0.1-10 keV; Parmar et al. 1997) and Medium- 
Energy Goncentrator Spectrometer (MEGS; 1.3-10 keV; 
Boella et al. 1997) on-board BeppoSAX are presented. The 
MEGS consists of three identical grazing incidence tele- 
scopes with imaging gas scintillation proportional coun- 
ters in their focal planes. The LEGS uses an identical con- 
centrator system as the MEGS, but utilizes an ultra-thin 
(1.25 fiTo) entrance window and a driftless configuration 
to extend the low-energy response to 0.1 keV. The fields 
of view (FOV) of the LEGS and MEGS are circular with 
diameters of 37' and 56', respectively. In the overlapping 
energy range, the position resolution of both instruments 
is similar and corresponds to 90% encircled energy within 
a radius of 2'.5 at 1.5 keV. At lower energies, the encircled 
energy is proportional to E~°-^. The LEGS 0.1-10 keV 
background counting rate is 9.7 x 10^^ arcmin^^ a~^. 

The region of sky containing IE 2259+586 and 
G109. 1-1.0 was observed by BeppoSAX between 1996 
November 16 03:32 and November 17 05:19 UTG. Good 
data were selected from intervals when the minimum ele- 
vation angle above the Earth's limb was >4° and when the 
instrument configurations were nominal using the SAX- 
DAS 1.1.0 data analysis package. This gives exposures of 
53.4 ks for the MEGS and 15.2 ks for the LEGS, which was 
only operated during satellite night-time. The LEGS and 
MEGS images, smoothed using a Gaussian filter of width 
1'.5, are shown in Fig. |l|. IE 2259+586 is visible as a point 
source close to the center of both images. The position 
of IE 2259+586 derived by summing the LEGS and the 
MEGS data is 23'' 01"" 06!8, 58° 52' 16" (J2000.0) with 
a 68% confidence uncertainty radius of 30". This is con- 
sistent with the Einstein High Resolution Imager position 
of Fahlman et al. (1982). The SNR is visible as extended 
emission throughout both images. 

Separate spectra were extracted for the pulsar, SNR 
shell and X-ray lobe regions. All spectra were rebinned 
to have >20 counts in each bin to allow the use of 
statistics. Background subtraction was performed using 
standard blank field exposures. 

3. Results 

3.1. The IE 2259+586 spectrum 

LEGS and MEGS spectra were obtained centered on the 
pulsar position using an extraction radius of 4'. The 
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Fig. 2. The best-fit power-law and blackbody fit to the 
pulsar spectrum (top) and the best-fit power-law model 
fitted to the pulsar spectrum with 0.19 of the SNR shell 
spectrum subtracted (bottom). The smaller panels give 
the residuals in counts keV^"'. The models are de- 
scribed in the text and summarized in Table 

IE 2259+586 count rates above background are 0.28 and 
0.55 s~-^ in the LEGS and MEGS, respectively. Examina- 
tion of the LEGS spectrum reveals that the pulsar is only 
detected between 0.5 and 8.0 keV and data outside this 
range are excluded. Similarly, the MEGS fit is restricted 
to the energy range 1.65-10 keV. 

The combined LEGS and MEGS spectrum of the pul- 
sar was first fit using a power-law model with a ~ 
4.49 ± 0.03 and low-energy absorption of (1.46 ± 0.05) x 
10^^ atoms cm~^ yielding a oi 449 for 266 degrees of 
freedom (dof). The photo-electric absorption coefficients 
of Morisson & McGammon (1983) and the solar abun- 
dances of Anders & Grevesse (1989) were used for all fits. 
Examination of the residuals reveals significant structure 
below -1.0 keV. A better fit (x^ of 271 for 264 dof) is ob- 
tained when a blackbody component is added. The spec- 
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Fig. 1. LECS (left panel) 0.5-2.0 keV and MECS (right panel) 2.0-8.0 keV smoothed images of the region containing 
IE 2259+586 and G109. 1—1.0. The contours levels represent 1, 2, 3, 6, 12, 25 and 50% of the peak flux of each image. 
The "cut-outs" at the upper-left and lower-right of the MECS image are due to the removal of instrument calibration 
source events. The softness of the SNR and lobe can be seen by the greater asymmetry of the LECS contours, compared 
to the MECS. The extraction region used for the SNR shell spectrum is indicated by the heavy dotted line, minus the 
box shaped region (which is the extraction region for the lobe spectrum) 



tral parameters listed in Table |l| are similar to those ob- 
tained using the same spectral model with ASCA spectra 
by Corbet et al. (1995) and with a combination of ASCA, 
ROSAT, and BBXRT spectra, by Rho & Petre (1997). The 
remaining residuals near 1 keV may indicate an even more 
complex spectral shape. 

In order to investigate an alternative explanation for 
the low energy residuals, the spectrum of the SNR shell 
multiplied by a scale factor, F, was subtracted from 
the IE 2259+586 LECS spectrum (the SNR hardly con- 
tributes in the MECS energy range, see Fig. The fit 
with a single power-law was then repeated and the best- 
fit value of F of 0.19 ± 0.01 determined by minimizing x^. 
This gives a of 296 for 265 dof, and the amplitude of the 
residuals ^1 keV is similarly reduced as with the power- 
law and blackbody fit (Fig. H). The best-fit value of a is 
4.86±0.04 and Nh is (2.18±0.07) x 10^^ ^toms cm'^. The 
ratio of the areas of the pulsar and SNR shell extraction 
regions (39 and 361 arcmin^, respectively) is 0.11. This is 
comparable with the value obtained for F of 0.19, given the 
observed variations in surface brightness of the SNR, the 
presence of the X-ray lobe, and the extended (2-3' radius) 
emission around IE 2259+586 reported in Rho & Petre 
(1997). Although the fit quality is not as good as with the 
power-law and blackbody model, we cannot exclude the 
possibility that some, or all, of the residuals below ^1 keV 



are caused by the contribution of the SNR that lies within 
the pulsar's extraction region. The IE 2259+586 spectral 
fit results are summarized in Table |l|. The 2-10 keV lumi- 
nosity of IE 2259+586 is 3.3 x lO^** erg s'^ for a distance 
of 4 kpc, identical to the value in Corbet et al. (1995). 

3.2. IE 2259+586 pulse timing and phase resolved 
.spectroscopy 

The MECS counts were used to determine the 
IE 2259+586 pulse period, after correction of their arrival 
times to the solar system barycenter. The data were di- 
vided into 14 time intervals (each with ^^2000 counts) and 
for each interval the relative phase of the pulsations de- 
termined. This was performed by folding the counts at 
half the pulse period value, in order to obtain light curves 
with a stronger modulation. The phases of the 14 time 
intervals were then fitted with a linear function giving a 
best-fit period of 6.978914 + 0.000006 s. The IE 2259+586 
light curve (Fig. H) shows a double-peaked profile with the 
amplitude of the second peak about half that of the main 
peak. 

A set of four phase-resolved spectra of the pulsar were 
accumulated, approximately coinciding with the peaks 
and valleys of the pulse profile. These spectra were fit 
with the same power-law plus blackbody model as used 
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Table 1. Spectral fit results for IE 2259+586. Uncertain- 
tics arc given at 68% confidence 



Parameter 


Value 


Pulsar (power-law model) 




a 


4.49 ± 0.03 


Nh (10^^ atoms cm~^) 


1.46 + 0.05 




449/266 



Pulsar (power-law & blackbody model) 
a 3.93 + 0.09 

kTbb (keV) 0.44 + 0.01 

Blackbody radius" (km) 3.3+E!:i 
Nh (10^^ atoms cm"^) 8.7 + 0.5 

271/264 

Pulsar-0.19.SNR (power-law model) 
a 4.86 + 0.04 

Nh (10^^ atoms cm"^) 2.18 + 0.07 

296/265 

"For a distance of 4 kpc 



in Sect, ^j, with Nh fixed at the phase- averaged best- 
fit value. There are insufficient counts to simultaneously 
constrain both the power-law and blackbody components. 
Initially, the blackbody spectral parameters were fixed at 
their phase-averaged best-fit values and only the power- 
law parameters allowed to vary. The fits were then re- 
peated with the power-law parameters fixed while the 
blackbody parameters were allowed to vary. With the lat- 
ter approach, the two fits in the valleys are unacceptable 
with x^'s of 2.9 and 2.3 for 117 dof. This is because the 
power-law component contributes too much flux, even if 
the contribution from the blackbody is set to zero. 

The best-fit values of a obtained with the first ap- 
proach are shown in Fig. |^ and reveal a small phase de- 
pendence. This variation corresponds to a of 9.4 for 
3 dof with respect to a constant value. Although this is 
significant at >99% confidence, we cannot be certain that 
there are real variations in a, since there is no clear corre- 
lation with the flux (as might be expected). In addition, 
small uncertainties in background subtraction, or a con- 
tribution from the SNR could cause such an effect, and so 
we prefer to set a limit of ±0.2 to any phase-dependent 
change in a. 

3.3. The G109. 1-1.0 spectrum 

In order to accumulate the spectrum of the SNR shell 
a complex shaped extraction region consisting of a circle, 
delimited by two straight lines with the pulsar and jet-like 
X-ray lobe cut out was used (see Fig. |l|). Since the spectra 
of both the shell and the lobe are soft (there is little flux 
^2 keV), only LEGS data were used. The effective area 
of the LEGS depends on source position within the FOV 
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Fig. 3. The pulse profile (top panel) of IE 2259+586 in 
the 0.5-10 keV energy range. The bottom panel shows 
the values of a. The values are repeated for clarity 

and the appropriate response matrix was determined by 
using 11 point sources whose position and relative inten- 
sities were chosen to mimic the observed count distribu- 
tion. Examination of the extracted spectrum shows that 
the SNR is only detected between 0.25-4.5 keV and data 
outside this range are excluded. 

The spectrum of the SNR shell was fit with the Non- 
Equilibrium Ionization (NEI) plasma emission model (plus 
absorption) implemented in V. 1.10 of the SPEX package. 
Fits were performed with freely varying O, Ne, Mg, Si, S, 
Fe and Ni abundances. A one-component NEI spectrum 
satisfactorily describes the spectrum with a of 65 for 66 
dof (see Table ||). The spectrum of the X-ray lobe was sep- 
arately extracted and analyzed, using the region indicated 
in Fig. ||. Since the size of the extraction region is compa- 
rable to the standard LEGS extraction region, no correc- 
tion for source extent was applied to the response matrix. 
A background spectrum was obtained from the standard 
blank field exposures using the same extraction region as 
for the source spectrum. When the X-ray lobe spectrum 
is fit with the same NEI model, the best-fit values are in 
all cases within Icr of those of the shell spectrum. This 
is in agreement with the results of Rho & Petre (1997) 
and supports the view that the two regions are physically 
related. 

4. Discussion 

The region of sky containing the IE 2259+586 pulsar and 
the G109.1-1.0 SNR is complex. The 1E2259+586 spec- 
trum is best described by the sum of a power-law and a 
blackbody, confirming the results of Gorbet et al. (1995). 
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Table 2. Spectral parameters for the one-component NEI 
fit to the SNR shell spectrum. Uncertainties are given at 
90% confidence (Ax^ = 2.71). A distance of 4 kpc is as- 
sumed. The X-ray lobe has consistent spectral parameters, 
except for an emission measure (uc uh V) of 27 x 10^® cm~'^ 



Parameter 




Best-fit 
value 


90% conf. 
range 


Nh (10^^ atoms cm 




6.9 


5.7-7.6 


kT (kcV) 




0.95 


0.68-1.60 


ncHH V (lO'^^cm-^ 


] 


63 


27-120 


Hot (10^cm"^yr) 




1.2 


0.7-2.4 


[O/H] 




2.0 


0.5-5.0 


[Ne/H] 




1.1 


0.0-1.2 


[Mg/H] 




2.4 


1.3-3.9 


[Si/H] 




1.1 


0.7-1.8 


[S/H] 




0.8 


0.1-2.0 


[Fe/H] 




1.9 


1.3-2.9 




G109. 
and 


1 -1 .0 SNR spectrum 
best— fit NEI mode 
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Fig. 4. The G109. 1—1.0 shell spectrum, together with the 
best-fit single-component NEI model 



These authors suggest that at least part of the power- 
law component could originate from a synchrotron neb- 
ula which may be visible around IE 2259-1-586 in ROSAT 
images (Rho & Petre 1997). A similar spectral decompo- 
sition has also been reported for two other "anomalous 
pulsars", 4U 0142-^61 and IE 1048.1-5937 (White et al. 
1996; Corbet & Mihara 1997) and has been interpreted as 
evidence for quasi-spherical accretion onto isolated neu- 
tron stars after common envelope evolution and spiral-in 
(Ghosh et al. 1997). In this case the accretion flow has two 
components. A low-angular momentum component giving 
rise to the blackbody emission from a large fraction of the 
neutron star surface, and a high-angular momentum one 
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Fig. 5. Pulse period history of IE 2259-1-586. See Baykal 
& Swank (1996) for the measured values, except for the 
BeppoSAX result reported here 



forming an accretion disk responsible for the power-law 
emission. 

We caution that an alternative explanation for the ob- 
served two-component spectrum cannot be excluded for 
IE 2259-1-586. If a plausible fraction of the emission arises 
from the part of the SNR within the pulsar's extraction 
region, then an acceptable fit can be obtained with a sim- 
ple power-law model. This implies that the derived black- 
body radius for IE 2259+586 should be regarded as an 
upper limit. We note that with the power-law fit, the pul- 
sar's Nh of >2.1 X 10^^ atoms cm~^ is significantly greater 
than obtained for G109. 1 — 1.0 of <7.6 x 10^^ atoms cm~^. 
This may indicate that the pulsar lies some distance be- 
yond the SNR and that the two objects are unrelated, 
or it may indicate the presence of absorbing material lo- 
cal to the pulsar. The Nh to the SNR is consistent with 
the galactic column in the direction of IE 2259+586 of 
7.4 X 10^^ atoms cm ^ (Dickey & Lockman 1990). 

The derived pulse period is consistent with the ex- 
trapolation of the long-term spin-down measured by Cor- 
bet et al. (1995) and Baykal & Swank (1996) as can be 
seen in Fig. |^. No large change in spin-down rate, as ob- 
served for example from 4U 1626—67 (Chakrabarty et al. 
1997), is evident. This implies that the accretion torque 
has remained approximately constant over at least the last 
19 years. 

In agreement with the results of Rho & Petre (1997), 
the LEGS spectra of the SNR shell and of the jet- like X-ray 
lobe are indistinguishable, supporting a common origin for 
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these two regions. We therefore refer to the summed spec- 
trum as "the SNR" . For an assumed distance of 4 kpc, 
the 15' radius of the SNR shell corresponds to 17.5 pc, 
or an emitting vohime of / x 6.6 x 10"''^ cm"*, where / is 
the filhng factor of the emitting plasma within a sphere. 
The mass of the SNR is then 60 x f°-^ Mq. The canon- 
ical filling factor for a strong shock is 0.25, although the 
absence of half of the remnant reduces this by a factor 
of 2. The dumpiness of the X-ray emission (with most 
of the emission coming from a few bright spots) implies 
an even smaller filling factor, and we adopt a value of 
0.1. The total emitting remnant mass is therefore ^15- 
20 M0. The remnant age is estimated to be 13,000 yr 
using the hydrodynamical model of Wang et al. (1992). 
However this estimate uses a remnant plasma tempera- 
ture of 0.4 keV, significantly lower than determined here. 
The higher temperature implies a faster expansion speed 
(900kms~^ rather than 590kms~^), and thus that the 
remnant is younger than previously estimated. The Wang 
et al. (1992) estimate is based on detailed numerical sim- 
ulations, and it is not straightforward to determine the 
age corresponding to the updated temperature. Similarly, 
Hughes et al. (1981) determine an age of 17,000 yr assum- 
ing an X-ray temperature of 0.17 keV. Given their analyt- 
ical approach, is possible to scale their estimated age to 
the current temperature determination to give a value of 
3000 yr. 

The low value derived for the emitting mass is consis- 
tent with a young remnant, which has not yet swept up 
large quantities of circumstellar material. The abundances 
derived from the X-ray spectrum are slightly higher than 
cosmic values, indicative of circumstellar material which 
has been mildly enriched either by the stellar wind of the 
progenitor in its late evolutionary stages, or by some mod- 
erate mixing-in of ejccta material. The value of the ioniza- 
tion parameter (1.2 x 10"^ cm^"^ yr) is indicative of strong 
non-equilibrium conditions, with an implied ionization age 
of 3000 yr, in good agreement with the scaled mmierical 
simulations of Hughes et al. (1981). The lack of strongly 
enriched material implies that the ejecta is not being di- 
rectly observed, and is also consistent with the lack of a 
strong wind of enriched gas originating from the progen- 
itor in its late stages. This is suggestive of a low-mass 
progenitor and a Type lb supernova. 

The LEGS fit results to the whole of the SNR are simi- 
lar to those of Rho & Petre (1997), who analyzed BBXRT 
and ROSAT spectra of a small part of the remnant, just 
south of the pulsar. Their fit with a NEI model in non- 
equipartition also shows strong non-equilibrium condi- 
tions, with a comparable (considered the differences in the 
instrumental responses and plasma emission models) age 
of 6700 yr, near-cosmic or slightly enriched abundances 
and an emitting mass of 85 Mq. 
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